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Intermediate-mass stars end their lives by ejecting the bulk of their enve- 
lope via a slow dense wind back into the interstellar medium, to form the next 
generation of stars and planets'^^. Stellar pulsations are thought to elevate 
gas to an altitude cool enough for the condensation of dust^, which is then 
accelerated by radiation pressure from starlight, entraining the gas and driv- 
ing the wind'SEE. However accounting for the mass loss has been a problem 
due to the difficulty in observing tenuous gas and dust tens of milliarcseconds 
from the star, and there is accordingly no consensus on the way sufficient 
momentum is transferred from the starlight to the outflow. Here, we present 
spatially-resolved, multi-wavelength observations of circumstellar dust shells 
of three stars on the asymptotic giant branch of the HR diagram. When 
imaged in scattered light, dust shells were found at remarkably small radii 
(< 2 stellar radii) and with unexpectedly large grains (~300 nm radius). This 
proximity to the photosphere argues for dust species that are transparent 
to starlight and therefore resistant to sublimation by the intense radiation 
fleld. While transparency usually implies insufficient radiative pressure to 
drive a wind*^, the radiation fleld can accelerate these large grains via pho- 
ton scattering rather than absorption*^ - a plausible mass-loss mechanism for 
lower-amplitude pulsating stars. 

We observed W Hya, R Dor and R Leo using aperture-masked'^'^ polarimetric inter- 
ferometry (see Figure 1), along with dust-free stars to verify our detection methodology. 
Figure 1 shows the ratio of the horizontally to vertically polarised visibilities (Vhoriz/Kert), 
plotted as a function of baseline azimuth and length. The dust-free star 2 Cen, which has 
no polarised flux, shows a constant ratio Vhoriz/Krert = 1-0 within uncertainties. However 
the dust-enshrouded star W Hya shows a strong sinusoidal variation of Vhoriz/Kert with 
azimuth, as expected from a resolved circumstellar scattering shell. By taking advantage 
of spherical symmetry, the much simpler baseline-dependent observable plotted in Fig- 
ure 2 was produced. A model was then fitted to the data to determine the dust shell 
radius and amount of light scattered by the shell at each wavelength (see Supplementary 
Information). These results are summarised in Table [l] Figure 3 shows the model image 
of the star and shell as seen in orthogonal polarisations for W Hya at 1.24/im, from which 
the model visibilities were derived. 

Scattered-light dust shells around the three AGB stars observed were found close to 
the star, with radii <2 Rgtar- This is in contrast to earlier models which place the shell at 
many stellar radii'^', but is consistent with some recent model^ and with interferometric'^ 
and polarimetric^ measurements. Based on typical elemental abundances and spectral 
observations, the composition of AGB dust shells is expected to be dominated by sili- 
(.g^^ggEHSm^ in the form of olivine (Mg2a;Fe2(i-a;)Si04) and/or pyroxene (Mga.Fe(i_a;)Si03) 
where < x < 1. The temperature of a grain is determined by its opacity - i.e. how 
strongly it absorbs the surrounding radiation field. Non-grey models'^ show that silicate 
dust which contains iron absorbs the stellar flux strongly (having high opacities at wave- 
lengths ~ 1/im where the energy distribution peaks) and so can only condense at distances 
greater than ~ 5 stellar radii. These iron-rich species could be accelerated via absorption 
of stellar radiation, but they form too far from the star to provide an efficient mass-loss 
mechanism for low-amplitude pulsators (semi-regular variables)'^. Our detection of dust 
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much closer to the star is instead consistent with the presence of iron-free sihcates such 
as forsterite (Mg2Si04) and enstatite (MgSiOa) which are almost transparent at wave- 
lengths of ~ Ifim. Such grains do not heat to sublimation, despite the intense radiation 
close to the star, but the same transparency also prevents the momentum transfer from 
starlight required to drive a wind. For some stars, a possible solution to this dilemma 
arises when very large grains are considered. 

The degree of scattering by dust depends strongly on wavelength and on the size of the 
particles. By analysing our multi-wavelength measurements using Mie scattering theory, 
the effective grain size and number of grains were determined. As shown in Figure 4, an 
effective grain radius of ~300nm was found. For grains of this size, the scattering opac- 
ity becomes very large, well beyond that arising from Rayleigh scattering when smaller 
particles are assumed. In this regime, the contribution to radiative acceleration via scat- 
tering, rather than by absorption alone, must be considered. Models show that grains 
exceeding a certain critical scattering opacity can drive a wind at high Mg condensation 
and that, for a star of temperature 2700 K, this critical scattering opacity is only ex- 
ceeded in a narrow range of dust grain radii around 300 nm^. It should also be noted 
that a narrow range of grain sizes, of the order of ~500nm radius, is predicted based on a 
self-regulating feedback mechanism: grain growth effectively halts once the critical size is 
reached, since the dust is then accelerated outwards and gas densities quickly decrease!^. 
This is consistent with observations of grains in the inter-stellar medium, which are dom- 
inated by silicates'^ and have similar grain size^^^*^. Wind driving due to scattering by 
Mg-rich silicates is consistent with the finding that mass-loss in these stars depends on 
their metallicity'^. While this model encounters difficulties in the case of stars with ex- 
tremely extended atmospheres, such as R Leo (due to the mass of the stellar atmosphere 
at and above the dust-forming layers being too high to allow sufficient acceleration'^, it 
provides a plausible explanation for the mass loss of semi-regular pulsating stars such as 
R Dor. Our observations provide direct evidence for a population of dust grains capable 
of powering a scattering-driven wind. 

The last column of Table [l] gives the mass of the dust that contributes to the observed 
scattering signal, assuming the shell to be thin and the dust grains to be forsterite of a 
uniform size. If full Mg condensation and solar abundances are assumed then the gas- 
to-dust ratio is ~600, yielding total shell masses of ~ 6 x 10""^ Mq, ~ 2 x 10"'' Mq 
and ~ 1 X 10"'' Mq for W Hya, R Dor and R Leo, respectively. Since the pulsation 
periods of these stars are ~1 year and the mass loss rates are ~ 1 x 10"^ Mq /year 20 22 
this implies that for stars with less extended atmospheres a significant fraction of the 
observed shell is ejected each pulsation cycle, consistent with the observed dust being 
part of an outflow. In the extended-atmosphere case, where full Mg condensation is not 
observationally supported, a possible alternative dust species is corundum (AI2O3), as 
discussed in the Supplementary Information. Further time-dependent grain growth and 
dynamical models will help elucidate the role played by light scattered from large-grained 
dust in the AGB star mass-loss process. 
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37.0 (fixed) 


0.022±0.004 






R Leo 


0.4 


1.04 


18.3±0.3 


29.9±0.4 


0.120±0.003 


~30D nmt 


~ 2 X 10-10 M0 



Table 1: Summary of fitted model parameters. The radii of the dust shells were 
found to be < 2i?star- The scattered fraction is the proportion of the total flux arising 
from scattering by the dust shell. The grain radius was obtained from fltting to multi- 
wavelength observations using Mie scattering (value for R Leo flxed at 300nm). The 
scattering-shell mass was calculated only for the observed dust (and not, for example, a 
distribution extending to small grains invisible to our technique). Stellar radii given are 
for a uniform disc. All three AGB stars - W Hya, R Dor and R Leo - were observed in 
March 2009, with additional observations of W Hya at L24yum and 2.06/im made in June 
2010. Although the photospheric diameter - and possibly the dust shell diameter - are 
expected to vary throughout the stellar pulsation cycle, the two sets of observations for 
W Hya were taken approximately one period apart, and so have been combined. These 
figures assume the dust to be an iron-free silicate (forsterite) grains of uniform size, and 
full Mg condensation. Full Mg condensation is a reasonable assumption for the stars 
with more compact atmospheres (e.g. R Dor) but is inconsistent with observed optical 
depths for stars with more extended atmospheres'^. In the event that there is also a 
population of small weakly scattering grains which do not show up in our data, these 
values represent lower limits with the total shell mass being higher. Furthermore, if the 
shell is geometrically extended then the true mass will be higher, as these calculations 
assume a thin shell. The uncertainties given are based on random errors and do not 
account for systematic errors such as those described here. Hipparcos parallaxes'^ and 
experimentally measured optical constants^- have been used. $ indicates visual phase, 
derived from the AAVSO International Database. ^Assumed grain size. 
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Figure 1: Polarimetric interferometry of W Hya at 1.24yum - Although hght scat- 
tered by each part of the circumstellar dust shell is strongly polarised, the integrated signal 
recovered with conventional polarimetry is zero for an unresolved spherically symmetric 
shell. In this study, aperture masking interferometry'SnQl (which converts the 8 m VLT 
pupil into a multi-element interferometer, using the NACOl^ instrument) allows access 
to the ~10 mas spatial scales required to resolve the shell, while polarimetric measure- 
ments (obtained by simultaneously measuring interferometric visibilities in orthogonal 
polarisations^^) allows light from the star and from the dust shell to be differentiated. 
Here, the ratio of the horizontally polarised visibilities (Vhoriz) to vertically polarised vis- 
ibilities (Kert) is plotted against baseline azimuth angle (corresponding to position angle 
on the sky). Colour encodes the baseline length (the longest being 7.3 m). \4oriz/Kert is 
a differential observable, which allows the cancellation of residual systematic errors and 
depends only on the fractional polarised scattered light signal. Panel a) shows the result 
for W Hya, an AGB star with a circumstellar shell; Vhoriz/Kert deviates from 1, varying 
sinusoidally with azimuth. This is the signal expected from a thin, spherically symmetric 
dust shell scattering the light from a central star. Different baselines exhibit varying 
amplitudes of this signal, encoding the spatial extent of the resolved structure. Data 
points have been repeated over 2 cycles. Error bars are la. The longest baselines (red) 
have poor signal-to-noise because they are close to the null where the visibility curve of 
the star is extremely low. Panel b) shows visibility data for the star 2 Cen, which has 
no circumstellar dust shell and hence no polarised signal from scattering; here, the ratio 
Vhoriz/Kert is ~1 for all azimuths. 
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Figure 2: Wavelength dependence of scattering for W Hya and R Dor - The 

azimuthaUy reduced Vhoriz/Kert for the AGB stars W Hya and R Dor, plotted against 
spatial frequency (baseline length/A), at multiple wavelengths. The functional form of the 
visibility ratio excursions exhibited by all three AGB stars are consistent, within uncer- 
tainties, with a simple spherically- symmetric shell. We were therefore able to significantly 
enhance our signal-to-noise in a quantitative analysis by reducing our two-dimensional 
visibility data to a one-dimensional function of the baseline length (corresponding to 
spatial frequency). This was achieved by dividing Vhoriz/Kert by the expected sinusoidal 
variation (characteristic of a spherical shell) at a fixed amplitude, resulting in the much 
simpler baseline-dependent observable plotted here. Points with la error bars are the 
observed data (binned) and the solid lines are the fitted model. The spatial frequency 
of the characteristic 'dip' varies with the radius of the dust shell: for a larger shell the 
dip occurs at lower spatial frequencies. The depth of the dip depends on the amount 
of scattered light, with a larger deviation from Vhoriz/Kert = 1-0 indicating that a larger 
fraction of the total fiux arises from light scattered by the shell. This is seen to decrease 
strongly at longer wavelengths as expected theoretically; the precise change in this quan- 
tity as a function of wavelength can be used to determine the dust grain radius using Mie 
scattering theory (see Figure 4). The fitted parameters for these quantities are included 
in Table [Tj The 2.06/im data have insufficient spatial resolution to constrain the shell 
size, so for models at these wavelengths the shell size was fixed to be consistent with the 
fitted size at shorter wavelengths. 
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Figure 3: Model image for W Hya with circumstellar shell viewed in horizon- 
tally and vertically polarised light - The white disc represents the uniform-disc star 
used in the model. A three-dimensional model of a star with a thin scattering shell was 
constructed, and the scattered intensity observed in each polarisation for each point on 
the shell was calculated using Mie scattering, finally yielding an image of the star and 
shell. Polarised visibilities were then derived from the model and fitted to the observed 
visibilities, to allow dust shell radius and scattered fraction to be determined. Details of 
the modelling process can be found in the supplementary information. Refer to Figure S2 
for a diagram illustrating how the polarised intensity distribution arises. 
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Figure 4: Grain size measurement - Grain-size fitted to the scattered light fraction as 
a function of wavelength, for W Hya. Error bars are Icr. The inset gives a blowup of the 
1.0-1.3 /im region. The solid line represents the fitted Mie scattering model (where grain 
size and grain number were fitted parameters), while the dashed line represents the best 
Rayleigh fit (grain size fixed to be below the Rayleigh limit). The data are inconsistent 
with Rayleigh scattering. The fitted Mie model yields an effective grain radius of 316±4 
nm. In reality a distribution of grain sizes may be present, for example a population of 
very much smaller particles would contribute only weakly and could lie undetected. Our 
data show that, regardless of the presence or absence of smaller grains, a population of 
large ~300nm grains is required. 



7 



Supplementary Information 



Modelling Details 

To determine the size of the dust shell and the proportion of the total flux arising from 
scattering by the dust shell, a model was created and fitted to the observed differential 
visibilities. The star was modelled as a uniform disk. While AGB stars may generally 
be better represented as a limb-darkened disk, it was found that using a limb-darkened 
disk model made a negligible difference to the fitted parameters describing the dust 
shell. Thus for the sake of simplicity, and to allow unambiguous reference to previously 
measured uniform-disk diameters, a uniform disk was used here. The assumption was 
imposed that the dust shell is of negligible thickness compared to its radius, and so 
was modelled as a thin spherically symmetric shell. While in reality the shell may be 
thick or form the inner edge of an outflow, due to the rapid fall-off with radius of both 
illumination and density (in the case of an outflow) the polarimetric signal is dominated 
by the bright inner surface, and hence our data is unable to constrain a more complex 
density profile. Furthermore, we are also unable to resolve the thickness of the shell. If 
the shell is significantly extended, then the true mass of the shell will be higher than the 
"scattering shell mass" estimated here. 

Next, images of the shell as seen in linearly polarised light were created. The following 
description refers to the angles labelled in Figure SI. The polarised intensity of light 
scattered by a dust grain towards an observer is a function of the angle between the 
ray incident upon the grain and the polarisation vector measured by the observer a, as 
well as wavelength, dust-grain size and refractive index of the material. Due to the close 
proximity of the dust to the star, the photosphere has a large angular size as seen from 
the point of view of a dust grain (represented by 6 and 0). Hence for each dust grain 
the polarised intensity of the scattered light seen by the observer was integrated over 6 
and 0, i.e. over the entire visible disc of the star. This scattered intensity was calculated 
using the University of Oxford Mie scattering code'^. It is assumed in this model that 
the shell is optically thin and multiple scattering does not occur. This calculation was 
repeated for each point in the dust shell in three dimensions, and then integrated along 
the line-of-sight of the observer to form a set of two-dimensional images of the shell. Once 
the image of the star was added to the image of the dust, a Fourier transform of each 
polarised image was used to produce the polarised visibilities of the star and shell, which 
could then be fitted to the observed data. 

Discussion 

In the extended-atmosphere case, where full Mg condensation is not supported by ob- 
servations, a possible alternative dust species to iron-free silicates is corundum (AI2O3), 
which is also largely transparent to the stellar flux. While the overall dust composition 
of AGB stars may be dominated by silicates, corundum may play an important role in 
dust formation and in radiative transfer processes-—. Gas temperatures at the observed 
dust shell radii derived from stellar models are low enough for both iron-free silicates and 
corundum to form (depending on the stellar phase)'^'^. Corundum generally forms ear- 
lier in the dust condensation sequence than silicates'^Sl and may form a condensation site 
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around which a sihcate mantle later condenses'^. If the dust observed here is corundum 
then this is consistent with other evidence for corundum existing as close as 1.5 stellar 
radii around AGB stars'^^^. The low abundance of Al means full condensation would 
yield an optical depth consistent with our data, and if the observed scattering shell is 
assumed to be entirely composed of corundum then a total shell mass ~20 times that 
given above is implied. In reality, the role of different dust species in the shell may be 
more complex. 




Figure SI: A diagram of the numerical model used to calculate the polarised intensities 
of scattered light. A description of the modelling procedure is included in the text of the 
supplementary information. The angle (j) is not shown for clarity. 
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Figure S2: A diagram sliowing liow tlie polarisation of starlight by a spherical dust 
shell arises. Unpolarised light emerges from the star and is scattered by dust in the 
spherical shell in the direction out of the page. This scattered light is linearly polarised, 
such that its polarisation is in a direction perpendicular to the plane defined by the light 
ray incident on the dust grain and the ray going from the dust grain to the observer. In 
this way, if the shell were viewed through a horizontal polariser, light reaching the viewer 
from the shell at ^ = 90° will have its full intensity while light coming from the shell at 
^ = 0° will have zero intensity (compare to Figure 3). 
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